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ABSTRACT 


In response to the escalating global demand for energy spurred by industrialization, nations are increasingly turning to alterna- 
tive fuel sources. Among these alternatives, hydrogen stands out for its remarkable efficiency and environmental benefits. With 
its potential to significantly reduce fuel consumption and air pollution, particularly in the transportation sector, hydrogen has 
garnered widespread attention as a promising solution to our energy challenges. This study focuses on the development of a 
dynamic model for an internal combustion engine powered by hydrogen (H2-ICE). By intricately dissecting and modeling the 
various components of the engine and their interactions, we aim to create a comprehensive simulation platform. This platform 
will enable us to accurately predict and analyze the performance of H2-ICE under different operating conditions. Through the 
utilization of advanced simulation tools like Matlab Simulink, we can validate the efficacy and reliability of our model, providing 
valuable insights into the behaviour of hydrogen-powered engines. Furthermore, the implications of this research extend beyond 
mere simulation. The developed model opens up avenues for further exploration and innovation in the field of hydrogen propul- 
sion. For instance, it lays the groundwork for the development of sophisticated control systems tailored specifically for H2-ICE 
applications. Additionally, the effectiveness of the proposed model in terms of air intake dynamics, pressure, temperature, fuel 
injection response, combustion efficiency, and overall engine performance, such as power output, torque, and engine velocity are 
illustrated through simulation results using Matlab Simulink. 
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1. Introduction 


With the increase in the global population, energy demand is ris- 
ing day by day. This not only depletes current energy reserves but 
also increases energy losses. Due to the limited potential of oil re- 
serves, it is necessary to increase the use of alternative fuels. The use 
of hydrogen in combination with petroleum-derived fuels in internal 
combustion engines can reduce harmful exhaust emissions from fos- 
sil fuels. That is why researchers continue to work on an economical 
and safe fuel that does not harm the environment [1 - 5]. In addition 
to hydrogen-enriched petroleum fuels in internal combustion en- 
gines, various companies continue to work on fuel cells. One of the 
most important properties of hydrogen is that it does not harm the 
environment at the end of its combustion. With the development of 
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technology, the importance of hydrogen will increase for reasons 
such as increased production, cost reduction, and portability. Hydro- 
gen is one of the important players in the clean energy sector. Hy- 
drogen is a portable, storable, usable, and sustainable energy carrier 
in many different areas. 

Over the past 30 years, extensive research has been conducted on 
the use of hydrogen as a fuel [6 - 8]. It is a directly usable energy 
source in the automotive sector with fuel cells. One of the most sig- 
nificant advantages of using hydrogen as a fuel in engines is that it 
reduces air pollution. Since there is no carbon in hydrogen fuel, there 
will be no emissions such as CO, CO», and HC at the end of com- 
bustion. When hydrogen energy is used in the transportation sector, 
which contributes to greenhouse effect, CO2 emissions will be re- 
duced. In a study conducted in Europe, it is estimated that by 2040, 
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new hydrogen vehicles will reduce CO2 emissions by 35%, while 
hydrogen-powered public transport vehicles will reduce CO2 emis- 
sions by 40%, with an average reduction in CO2 emissions of 44.8 
g/km [9]. 

Internal combustion engines are the most significant factors lead- 
ing to an increase in environmental pollution [10 - 15]. The role of 
the internal combustion engine is to convert chemical energy into 
mechanical energy. In this engine, the system is renewed at each cy- 
cle. The system is in contact with only one heat source (the atmos- 
phere and the fuel). Heat is produced by combustion in a variable 
volume chamber. It serves to increase the pressure within a gas fill- 
ing this chamber (this gas is also initially composed of fuel and oxi- 
dizer: air). This increase in pressure results in a force exerted on a 
piston, a force that transforms the translational motion of the piston 
into the rotational motion of the shaft (crankshaft). The operation 
takes place in 4 stages: intake; compression; ignition (combustion); 
exhaust [11]. 

The most popular internal combustion engines are classified into 
two categories according to the ignition technique of the air-fuel 
mixture: spark ignition engines (gasoline engines) and compression 
ignition engines (Diesel engines). 

A Diesel engine operates on the principle of self-ignition of the 
fuel. The richness of the mixture characterizes the engine load. The 
higher the load, the greater the quantity of injected fuel, while the 
airflow remains constant. When the piston rises in the combustion 
chamber (compression phase), it greatly increases the temperature 
and pressure, allowing the air-fuel mixture to self-ignite (the self- 
ignition potential of a fuel is characterized by the cetane number). 

In a gasoline engine, the load is characterized by the airflow that 
varies depending on the load control. For an optimal compromise 
between performance, consumption, and emissions, gasoline en- 
gines operate around stoichiometry. Moreover, this promotes effi- 
cient operation of the three-way catalytic converter (for the three 
main pollutants: CO, NOx, and HC). Indeed, the efficiency of this 
type of catalyst significantly degrades when deviating from stoichi- 
ometry. 

Research indicates that most of the CO, HC, and NOx compounds 
produced as a result of the combustion of fossil fuels come from gas- 
oline and diesel engines. Harmful waste such as sulfur dioxide, lead, 
and soot are also present in the exhaust emissions of internal com- 
bustion engines [11], [16]. 

The sustainable use and production of hydrogen, which does not 
harm the environment, will be one possible solution. A hydrogen in- 
ternal combustion engine vehicle is a type of hydrogen vehicle using 
an internal combustion engine. Hydrogen internal combustion en- 
gine vehicles are different from hydrogen fuel cell vehicles (which 
use hydrogen through electrochemical means rather than combus- 
tion). Instead, the hydrogen internal combustion engine is simply a 
modified version of the traditional gasoline internal combustion en- 
gine. Zero carbon means that no CO is produced, thus eliminating 
the primary greenhouse gas emission from a conventional petroleum 
engine. 

The characterization of an engine significantly affects its perfor- 
mance, leading to numerous studies aimed at optimizing engine per- 
formance. For instance, in [17], the surface of the combustion cham- 
ber is analyzed to understand heat transfer correlations. Similarly, 
the piston crown in relation to various combustion mode strategies 
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is investigated in [18]. The influence of fuel system variations on 
performance and emission characteristics is explored in [19], while 
[20] examines how combustion chamber geometry and fuel supply 
system variations can economize fuel consumption and reduce ex- 
haust emissions. The impact of blended fuels with split injections on 
combustion and emission characteristics is analyzed in [21]. Addi- 
tionally, [22] studies the effect of piston crown shape and the posi- 
tioning of spark plugs and fuel injectors on fuel system control and 
emissions. Research in [23] delves into the use of hydrogen as a fuel 
for internal combustion engines, presenting a detailed analysis. Ex- 
perimental studies focused on internal combustion engines fueled 
with hydrogen are also conducted in [24 - 25], offering valuable in- 
sights into the practical application and benefits of hydrogen as an 
alternative fuel. These comprehensive investigations collectively 
contribute to the advancement of engine design and fuel system op- 
timization, highlighting the potential for improved efficiency and re- 
duced environmental impact. 

Motivated by the aforementioned studies, this research delves into 
the utilization of hydrogen as a fuel in an ICE, exploring a promising 
avenue towards reducing emissions and advancing sustainability in 
automotive propulsion systems. By leveraging the dynamic model 
initially developed for gasoline-powered ICEs, this study aims to 
adapt and extend it to accommodate the unique characteristics and 
operational parameters associated with hydrogen combustion. The 
goal is to construct a comprehensive and accurate dynamic model of 
an H2-ICE, capable of simulating its performance under various op- 
erating conditions. 

The proposed model facilitates further investigation into the ef- 
fects of different components, including engine characterization, 
fuel characterization, and various methods of mixing fuel. For in- 
stance, by examining how hydrogen fuel impacts the combustion 
process, insights can be gained into optimizing fuel injection timing, 
air-fuel ratio, and ignition parameters to achieve efficient combus- 
tion and minimal emissions. Additionally, the model allows for the 
analysis of different hydrogen blending techniques and their effects 
on engine performance and emissions. 

This study's dynamic model also provides a platform to explore 
the thermal management challenges unique to hydrogen combustion, 
such as its higher combustion temperature and faster flame speed 
compared to gasoline. Understanding these dynamics is crucial for 
designing effective cooling systems and ensuring engine durability. 
Moreover, the model can simulate the impact of hydrogen's lower 
energy density on fuel storage and delivery systems, offering in- 
sights into optimizing fuel tank design and refueling infrastructure 
for hydrogen-powered vehicles. 

The effectiveness and fidelity of our proposed model will be rig- 
orously assessed and validated using Matlab Simulink, a versatile 
platform widely employed for dynamic system modeling and simu- 
lation. Through meticulous calibration and validation against empir- 
ical data, we seek to ensure that our model accurately captures the 
dynamic behavior and performance nuances specific to hydrogen- 
fueled ICEs. This validation process is essential for establishing the 
reliability and predictive capability of our model, thus instilling con- 
fidence in its utility for subsequent analyses and optimization efforts. 


2. Dynamic model of an H2-ICE 


Four-stroke internal combustion engines are the preferred choice 
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for conventional and hybrid vehicles due to their inherent fuel effi- 
ciency advantages over their two-stroke counterparts. This effi- 
ciency arises from the distinctive four-stroke engine cycle, where a 
power stroke occurs once every two revolutions of the crankshaft. 
This stands in contrast to two-stroke engines, which feature more 
frequent power strokes happening every revolution, resulting in 
higher fuel consumption per unit of power produced. Consequently, 
four-stroke engines boast reduced fuel consumption, emissions, and 
torque output compared to their two-stroke counterparts. 

Key components of the four-stroke engine include the intake and 
exhaust manifold, intake and exhaust valves, spark plug, piston, 
coolant, cylinder, crankcase, connecting rod, and crankshaft. Each 
of these components plays a pivotal role in facilitating the engine's 
operation and performance. This study centers on the development 
of a comprehensive model for the four-stroke internal combustion 
engine, leveraging foundational research outlined in prior studies. 
Figure | provides a graphical representation of the fundamental op- 
erating principle of an internal combustion engine. It illustrates the 
sequential processes of intake, compression, power generation 
within the combustion chamber, and the subsequent transmission of 
torque to the crankshaft to drive the engine. Through meticulous 
modeling and analysis, this research aims to enhance our under- 
standing of the intricate dynamics inherent to four-stroke engines, 
ultimately optimizing their performance in various automotive ap- 
plications. 
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Figure 1. Representative operating diagram of the SI engine. 
2.1. Intake Air System 


The Intake Air System in an ICE is a crucial component respon- 
sible for delivering the appropriate amount of air to the engine’s 
combustion chambers. Proper air intake is essential for efficient 
combustion and optimal engine performance. This system represents 
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the flow limit, which typically includes ambient temperature, pres- 
sure, specific enthalpy, and mass fraction (in the form of dry air). 
The ambient temperature and pressure (ambient temperature) can be 
constant or provided externally. The specific enthalpy is calculated 
as follows: 

(1) 


Nair = Cogir Tat 


Where c,_ i, is the heat coefficient at constant and T,;, is the am- 
bient temperature equal to 300K. 

The mass fraction of dry air represents the percentage of different 
components in the inlet air, which are: Oxygen (O2), Nitrogen (Nz), 
Unburned fuel (this component only appears after combustion), Car- 
bon dioxide (CO2), Water vapor (H2O), Carbon monoxide (CO), Ni- 
tric oxide (NO), Nitrogen dioxide (NO2), Nitric oxide and nitrogen 
dioxide (NOx), Particulate matter, Air, Exhaust gases. Each compo- 
nent's mass fraction contributes to the overall composition of the in- 
let air, playing a crucial role in determining the system's behavior 
and performance characteristics. 


2.2. Throttle Body 


The throttle body is the part of the air intake system that controls 
the amount of air flowing into the engine, in response to driver ac- 
celerator pedal input in the main. The throttle body is usually located 
between the air filter box and the intake manifold, and it is usually 
attached to, or near, the mass airflow sensor. The mass flow rate 
through the throttle is determined as follow [16][26]: 

The gas mass flow through the orifice: 


_ Acre Pupstr w(P : ) 
_ ratio 
Tupstr 


Thorp (2) 


Where P,a:i.is the ratio of downstream pressure (Paowser) and up- 
stream pressure (Pypstr)- 


Y(Pratio) 
y+1 4 
( 2 ie ifs P 4 ( 2 yr 
are Uf: Fratio yt+1 
2 2 Jt. (3) 
Y Po. Y y-1 
= y-1 > if (- a :) < Pratio < Piim 
ratio 
Pratio — 
Pim = 


Where Pypser = Pape = 101325 Pa is the ambient pressure. 
Y(P_atio) is the based follow correlation and is the function of 
P,atio Calculates the based on the different conditions of flow. 


(4) 


P _ Paowstr 
ratio — P 
upstr 


Where Pi; is the pressure ratio limitation, y = 1.4is the ratio of 


specific heats. The air mass flow rate (77g;,) then can be calculated 
as follows: 


(5) 


Mair = MorfVupstr air 


Where Yunstr_air 18 the upstream species mass fraction. 
Ae ffs the effect area, and is calculated as follows: 
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T 
Aer f = q ene” Cenr cnr) (6) 
90 (7) 
Otnr = Oce_tnr 100 


Where 6;,, is the open angle of the throttle (in degree), 
Oct thr is the percentage of throttle body that open, Dyn, is the 
throttle diameter at opening, C;;(@:n7-)is the discharge coefficient. 

The heat flow rate qo;¢ of is calculated as follows: 

dorf = Morf hupser (8) 
Where hupstris the upstream specific enthalpy (hgi,). 
2.3. Intake Manifold 


To calculate the fuel flow in the hydrogen (H2) internal combustion 
engine, one utilizes the characteristics of the fuel injector along with 
the fuel injector pulse-width. The model for fuel flow is expressed 
as: 
ia = NSinjPwingNeyt (9) 
fuel“ “¢,,1000 * 60 
Where N represents engine speed in rpm, Sjnj; denotes the fuel 
injector slope, Ny, signifies the number of engine cylinders, 
Cys Stands for the crankshaft revolutions per power stroke. 
The fuel heat flow rate is subsequently determined by: 
dfuel = Muetupstr (10) 
In contrast to gasoline spark-ignition (SI) engines, where gasoline 
is the sole fuel, in the case of hydrogen (H2) engines, the fuel also 
serves as the intake gas. Therefore, the intake gas mass flow rate and 
its associated fuel heat flow must be corrected. The corrected intake 
gas mass flow rate (Mn intx) is calculated as the sum of the fuel 
mass flow rate (Myer) and the gas mass flow rate (gas). 
Min intk = Mpuel + Mgas (11) 
Similarly, the corrected intake gas heat flow (in intx) is determined 
by summing the fuel heat flow (qfye,) and the gas heat flow (orf). 


Win_intk = fuel + orf (12) 


Subsequently, the species mass fractions are reformulated as follows: 


y _ Mpuel (13) 
unbrn_fuel encanate 
™Mgas (14) 
Yair =. _ 
Min intk 
m 15 
Yor = 0.233-— ve 
Min _intk 
(16) 


m 
Veo. = 0.767 — 
Min intk 


In this study, the other mass fractions remain zero. These formu- 
lations allow for a comprehensive understanding of the fuel and gas 
flows within the hydrogen internal combustion engine, facilitating 
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accurate modeling and analysis of its performance characteristics. 
The intake manifold serves as a critical component of the air intake 
system, responsible for regulating the pressure and temperature of 
the air delivered to the cylinder. The changes in temperature and 
pressure within the intake manifold are dynamically determined by 
employing a constant volume chamber model containing an ideal 
gas. This analysis leverages the continuity equation and the first law 
of thermodynamics to accurately characterize the behavior of the in- 
take manifold. The rate changes in temperature and pressure are de- 
termined by implementing a constant volume chamber containing an 
ideal gas, the continuity equation and the first law of thermodynam- 
ics is used in this case as follows: 


ATintake R Tintake ((4 (17) 
= in_intk 
dt Cp_airVenPintake 
™ ‘intake Cytlten siti) 
~ (Gout inet 
™ ‘intake Gptcucsaak) ) 
dPintake ia Pintake AT intake (18) 
dt Tintake dt 
RT intake : : 
+ i, (Thin intk — Mout ee) 
ch 7 ~ 


Where R is Ideal gas constant, V,, is Chamber volume, 
Tintake @0d Pintake are the chamber absolute temperature and 
pressure respectively, Gout intx 1S the output heat flow given by Eq. 
17, and Mout intz 1S the mass flow in the chamber given by Eq. 25. 
Hence the specific enthalpy of the intake manifold is calculated as: 
Rintake = Cp air intake (19) 
Using the conservation of mass for each gas constituent, the mass 
fraction of gas in the chamber can be determined as follows: 
AYintake _ RTintake 


dt Von Pintake 


(20) 


((vin inti Titin inek 
— Yintake Mout intic) 


+ Yintake (Coe = fitgaeseh)) 


These equations provide a comprehensive framework for under- 
standing the dynamic behavior of temperature, pressure, and gas 
composition within the intake manifold, facilitating accurate model- 
ing and analysis for optimizing engine performance. 

To determine the output heat flow (dout intk), We employ a heat 
transfer model depicted in Figure 2. This model accounts for both 
convective and conductive heat transfer processes occurring within 
the intake manifold. The output heat flow is calculated as the sum of 
convective (Q; cony) and conductive (Q; cong) heat transfer compo- 
nents. 

(21) 


out_intk = Qi = Qi conv = Qi cond 


Engineering Perspective 4 (3): 108-118, 2024 


Q.T.Dam et al. 
Q, = Q> conv = Q2 cond (22) 
Q2 
Ql 


Figure 2. Heat transfer model 


Firstly, Q1,convtepresents the convective heat transfer from the in- 
take gas to the interior wall of the manifold. It is computed using the 
equation: 
Qi conv = Cg int iiss )Armecansl Pages ad lie) (23) 
Here, cq ine represents the convective heat transfer coefficient, 
Min intk denotes the mass flow rate of the intake gas, Aint cony 18 the 
convective surface area, Tj, gas 1S the temperature of the intake gas, 


and 7, in. represents the temperature of the interior wall of the in- 


take manifold. 

Secondly, Q1,cona denotes the conductive heat transfer from the in- 
terior wall of the manifold to the mass of the wall itself. It is given 
by: 

KintAint cona (24) 


Qi cond oa (Todas J Tes) 


Dint cond 


Where kj; is the thermal conductivity of the wall material, 

Aint cona tepresents the conductive surface area, Dint cong is the 
thickness of the wall, and Tis; denotes the temperature of the wall 
material. 
Similarly, the heat transfer from the exterior wall of the manifold to 
the external environment is represented by Q2cony and Q2 cond: 
These components account for convective and conductive heat trans- 
fer processes, respectively, and are calculated using analogous equa- 
tions, and are given as follows: 


Q2,conv = Cq ext (ExtnIFlwVel) Aext conv (Tw_ext (25) 
~ Txegas) 
RextAext cond (26) 


= Toe) 


Q2 cond = (lies 


Dext cond 


Finally, the rate of change in temperature of the wall material 
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(Fa85) is determined by the difference between the input and output 


heat flows divided by the thermal capacity of the wall material as 


follows: 
AT mass _ 


dt 


Qi — Q2 


Cowall Mwall 


(27) 


By incorporating Eq. (13) to Eq. (19) into this heat transfer model, 
we can accurately derive the output heat flow (qoute intx), providing 
insights into the thermal dynamics of the intake manifold and facili- 
tating optimization of heat management strategies for enhanced en- 
gine performance. 

The Engine Cooling System serves several critical functions 
within the internal combustion engine (ICE) setup. Primarily, it acts 
to dissipate excess heat generated during engine operation, thereby 
preventing overheating and potential damage to engine components. 
Additionally, the cooling system plays a pivotal role in maintaining 
the engine at its optimal operating temperature, ensuring peak effi- 
ciency and performance. Furthermore, during engine startup, the 
cooling system facilitates the rapid attainment of the desired operat- 
ing temperature, promoting smoother and more efficient engine op- 
eration from the outset. 

In this study, our focus is not on designing the Engine Cooling 
System per se, but rather on integrating it seamlessly into our ICE 
system to form a closed-loop system. This integration is crucial for 
ensuring the overall functionality and efficiency of the ICE setup. 
The temperature of the engine after undergoing the cooling process 
is governed by the following equation: 

AT cooling 
dt 


(28) 


_ Mout intk (Gp: te Vitae O a nip oxi) + LHVMpye1 — TH 


Ceng 


Where LHV is the fuel lower heating value, w is the engine 
speed in rad/s. Ceng = 40000 is the heat capacity. 

The H2-ICE model employed in this study is based on the Spark- 

Ignition (SI) engine model, with hydrogen (Hz) serving as a substi- 
tute for gasoline. This choice reflects a shift towards cleaner and 
more sustainable fuel alternatives, aligning with the broader objec- 
tive of reducing emissions and promoting environmental sustaina- 
bility in automotive applications. 
By incorporating the Engine Cooling System into our H2-ICE model, 
we can ensure the holistic simulation and analysis of the entire en- 
gine system, enabling comprehensive assessments of performance, 
efficiency, and thermal management strategies. This integrated ap- 
proach lays the foundation for advancing the development and opti- 
mization of hydrogen-powered internal combustion engines for fu- 
ture automotive applications. 


2.4. Air Mass Flow Model 


To calculate the air mass flow in the engine chamber, one shall use 
a speed-density air mass flow model as depicted in Figure 3. The 
speed-density model uses the speed-density equation to calculate the 
engine air mass flow. The equation relates the engine air mass flow 
to the intake manifold gas pressure, intake manifold gas temperature, 
and engine speed. Consider using this air mass flow model in simple 
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conventional engine designs, where variable valvetrain technologies 
are not in use. 


Manifold 
Gas 
Pressure 
(MAP) 


Manifold Gas 
Temperature 
(MAT) ity = "Vt *Rair* 
Speed density = (MAP*Vd*N)/(Cps*Rair*MAT) 


Figure 3. Engine Chamber Speed-Density Air Mass Flow Model 


To determine the air mass flow, the speed-density air mass flow 
model applies these speed-density equations at the intake manifold 
gas pressure and gas temperature states. 

The mass flow (Mintaxe) at the intake port can be calculated as fol- 
lows [16]-[26]: 
_ = Pintake VaN n 
intake ah iene 60 /” 


(29) 


Where 7), is engine volumetric efficiency as given in Figure 4, di- 
mensionless, is a function of intake manifold absolute pressure 
(Pintaxe ) and engine speed (N) and calculated as follows: 
CosRTintaxe 60 


Pintake VaN 
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Figure 4. Engine volumetric efficiency 
Subsequently, the calculated air mass flow is multiplied by the mass 
fraction of air (Yintake air) in the intake mixture to determine the air 
mass (Mair intake): 


Mair intake = Vintake_air™intake 


(31) 
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Furthermore, the engine load (Leng), a crucial parameter indicative 
of the engine's operational state, is derived using the formula: 
Mair intake Cys * 6ORT sta 


VaN Psea 


(32) 
Leng = 


Where Tq, Pstq are the standard temperature and pressure re- 
spectively. 
The Air-Fuel ratio is calculated as: 
(33) 


Mair intake 


ARF = 


Meuel 


2.5. Torque Model 


In the SI engine the indicated torque means that the torque or power 
of the engine is evaluated in the scope of thermodynamics (pressure 
and volume of cylinder), not including any mechanical losses in the 
whole power development and transmission process conceptually il- 
lustrates pressure variation in a cylinder along with crankshaft rota- 


tion angle. This torque is calculated as follows [11]: 
= Pindicated * 60 
Tindicated ~~ 7 


(34) 


Where Pingicatea iS the indicated power and can be calculated as: 


Pindicated = Mpuel * LHV Ninaicated (35) 


Where YNindicated 1S the engine indicated load efficiency as well 
as a function of engine speed (N) and indicated torque (Tingicated) aS 
depicted in Figure 5. 


0.35 
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Tindicated (Nm) 
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Figure 5. Indicated torque efficiency 


Not all the torque produced in the cylinder (indicated power and 
torque) are available on the crankshaft. The available torque on the 
crankshaft is the difference between the indicated torque and friction 
torque, respectively, expressed as: 

Tprake — Tindicated Tfriction (36) 
The basic engine accessories include a water pump, oil pump, fuel 
pump, and valve shaft. Other accessories may depend on applica- 
tions, such as cooling fan, alternator, power steering pump, and air 
conditioner compressor. In an engine test, the accessories that are 
driven by the engine crankshaft should be clearly specified. All 
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power and torque requirements for driving accessories and overcom- 
ing frictions are grouped together and called the friction torque and 
is given as: 

Tfriction = Nprictionf (Leng N) (37) 
where Friction 18 the friction temperature factor as demonstrated 
in Figure 7, f(Leng,N) is the function of engine load (L,,,) and en- 
gine speed (N) as depicted in Figure 6. 
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Figure 8. Detailed Simulink model for ICE 


To illustrate the performance of the developed model, an explica- 
tive simulation will be conducted using Matlab Simulink, as shown 
in Figure 8. The simulation will utilize input references for the throt- 
tle valve open angle, given in Figure 9, and the fuel injector pulse- 
width, given in Figure 10, to demonstrate the model's performance. 
The scenario involves first regulating the throttle valve's open angle 
according to Figure 9, and then injecting hydrogen fuel into the en- 
gine following the values in Figure 10. This simulation aims to show 
how the throttle valve configuration affects overall system behavior, 
examining the relationship between throttle position and fuel injec- 
tor pulse-width in controlling the air and fuel intake. Expected out- 
comes include insights into air intake dynamics, pressure, tempera- 
ture, fuel injection response, combustion efficiency, and overall en- 
gine performance, such as power output, torque, and engine velocity. 
By integrating these findings, a comprehensive understanding of the 
throttle valve's influence on system behavior is achieved, providing 
a foundational step towards optimizing system design and enhancing 
operational efficiency. This detailed analysis ultimately validates the 
proposed H2-ICE model and offers valuable insights for further re- 
search and development in hydrogen-fueled automotive technolo- 
gies. 
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Figure 9. Throttle angle 
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Figure 10. Fuel injector pulse-width 


The ambient temperature and pressure are set up as 300 °K and 
101325 Pa respectively, the Heat coefficient at constant Cy gir = 
1005 J/(kg*K). We consider the dry air contains only oxygen and 
nitrogen with its mass fraction are given as Yo2 = 0.233 and 
Yn2 = 0.767. 

The downstream pressure within the intake manifold is calculated 
by applying Eq. (18), and the outcomes are meticulously presented 
in Figure 11. Subsequently, this calculated downstream pressure 
serves as a crucial parameter in determining the pressure ratio, as 
stipulated by Eq. (3), where the upstream pressure aligns with the 
ambient pressure. With this pressure ratio established, as graphically 
depicted in Figure 12, an essential step ensues in discerning the base 
flow correlation governing the gas mass flow through the orifice (2). 
This correlation is meticulously illustrated and analyzed in detail, of- 
fering valuable insights into the dynamics of gas flow within the sys- 
tem. 
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insights crucial for optimizing system performance and efficiency. 
Within this simulation, hydrogen fuel injection into the intake man- 
& ifold is implemented, with the corresponding fuel mass flow rate de- 
3 picted in Figure 15. These parameters play a pivotal role in regulat- 
: : ing the combustion process within the engine, ensuring optimal fuel- 
0 50 100 ts0 aif mixture ratios and combustion efficiency. By accurately model- 
Time (Sec) ing and analyzing the fuel injection dynamics, insights into engine 
Figure 11. Downstream pressure in the intake manifold performance, fuel consumption, and emissions can be gleaned. This 
comprehensive examination of fuel injection characteristics enriches 
1 our understanding of engine behavior under varying operating con- 
ditions and facilitates the refinement of engine control strategies for 
enhanced performance and efficiency. 
Utilizing the thermal dynamic equations of Eq. (21) to Eq. (27), 
the heat flow from the internal gas to a specified wall depth, termed 
fs = - tsp 2S the output heat flow, can be accurately determined. This crucial 
Time (Sec) parameter, which characterizes the heat transfer process within the 
Figure 12. Pressure Ratio system, is meticulously presented in Figure 16. Additionally, the 
heat transfer rate, representing the rate at which heat is exchanged 
06 ‘a _ — \f é 
between the gas and the surrounding walls, is quantified based on 
the same equations and visually depicted in Figure 17. 
0.2 These analyses offer valuable insights into the thermal dynamics 
o+¥__—— | of the system, shedding light on heat distribution, dissipation, and 
0 50 100 150 Overall thermal performance. By comprehensively understanding 
Time (Sec) the heat transfer processes, engineers can optimize system design, 
Figure 13. Base flow correlation enhance thermal efficiency, and mitigate thermal-related issues, 
thereby ensuring reliable and efficient operation under diverse oper- 
ating conditions. 
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Figure 14. Gas mass flow through the orifice 
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Given the determined output heat flow rate, it becomes feasible to 
ascertain the intake temperature, also referred to as the downstream 
temperature, as outlined in Eq. (17). This temperature parameter 
plays a pivotal role in regulating engine performance and combus- 
tion efficiency. By accurately modeling and analyzing the heat trans- 
fer dynamics within the system, insights into the variation of intake 
temperature under different operating conditions can be gained. Fig- 
ure 18 visually presents these temperature variations, offering a com- 
prehensive understanding of how heat transfer processes influence 


Time (Sec) 


Figure 16. Output heat flow rate 


Figures 13 present a comprehensive visualization of this base flow 
correlation, enriching our understanding of the intricate interplay be- 
tween pressure differentials and gas mass flow rates as presented in 
Figure 14. This detailed analysis not only elucidates the fundamental 
principles governing gas flow behavior but also provides practical 


115 


Q.T.Dam et al. 


Engineering Perspective 4 (3): 108-118, 2024 


intake temperature. This analysis enables engineers to optimize en- 
gine operation, ensuring that intake temperatures remain within de- 
sired ranges to promote efficient combustion and mitigate potential 
thermal-related issues. 

The mass fraction in the downstream of the intake manifold is de- 
termined using Eq. (20), which accounts for various factors influenc- 
ing the composition of gases within the system. By applying this 
equation, the simulation yields valuable insights into the down- 
stream mass fraction, representing the proportion of different gas 
components present in the intake manifold. Figure 19 visually illus- 


trates the results of this simulation, offering a comprehensive depic- 
tion of how the mass fraction varies under different operating condi- 
tions. This analysis enhances our understanding of gas composition 
dynamics within the intake manifold, providing essential infor- 
mation for optimizing combustion processes, controlling emissions, 
and ensuring efficient engine operation. By closely examining the 
downstream mass fraction, engineers can make informed decisions 
regarding fuel-air mixture ratios and combustion strategies, ulti- 
mately optimizing engine performance and reducing environmental 
impact. 
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Figure 20 depicts the gas mass flow rate within the chamber, 
providing crucial insights into the dynamics of gas flow through the 
engine system. This parameter serves as a fundamental metric for 
evaluating engine performance and efficiency. Subsequently, utiliz- 
ing the gas mass flow rate data, the indicated torque can be deter- 
mined according to Eq. (35), as illustrated in Figure 22. This torque 
represents the theoretical output torque generated by the engine, re- 
flecting the mechanical power produced during the combustion pro- 
cess. Conversely, Figure 22 presents the friction torque, which ac- 
counts for mechanical losses within the engine system. 

By comparing the indicated torque with the friction torque, engi- 
neers can assess the overall efficiency and mechanical integrity of 
the engine. Finally, Figure 23 shows the engine speed, providing val- 
uable insights into the rotational speed of the engine's crankshaft un- 
der varying operating conditions. Together, these parameters offer a 
comprehensive understanding of engine performance, enabling en- 
gineers to optimize design parameters, refine control strategies, and 
enhance overall operational efficiency. 
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4. Conclusions 


In conclusion, this study has centered on the examination of hy- 
drogen’'s potential as an alternative fuel for proposed hydrogen-pow- 
ered internal combustion engine (H2-ICE). Through the adaptation 
of a dynamic model originally designed for gasoline, a comprehen- 
sive dynamic model tailored specifically for ICEs fueled by hydro- 
gen has been successfully developed. This approach underscores the 
increasing significance of hydrogen as a sustainable energy carrier, 
presenting a promising pathway towards cleaner transportation so- 
lutions. Throughout the simulation results, the effectiveness of the 
H-ICE model was demonstrated using the throttle valve open angle 
and fuel injector pulse width as references. The model accurately re- 
flected the dynamic behavior of the engine under varying conditions, 
showcasing its potential for optimizing performance and efficiency. 
The findings underscore the model's capability to simulate real- 
world scenarios, providing valuable insights for future research and 
development in hydrogen-fueled engines. This study highlights the 
viability of hydrogen as an alternative fuel, contributing to the ad- 
vancement of sustainable and environmentally friendly automotive 
technologies. The exploration into the intricate dynamics of hydro- 
gen-powered internal combustion engines contributes substantially 
to ongoing research and development efforts aimed at advancing 
cleaner and potentially more efficient automotive technologies. By 
delving into the operational intricacies of these engines, this study 
lays the groundwork for future optimization and design endeavors 
within the realm of hydrogen-powered vehicles. Furthermore, this 
work not only addresses the pressing need for sustainable energy so- 
lutions but also opens up new avenues for innovation and advance- 
ment in the field of hydrogen vehicles. Through the successful adap- 
tation and validation of the dynamic model for H2-ICE, valuable in- 
sights into their operational characteristics, efficiency are provided, 
thus facilitating informed decision-making and strategic planning in 
the automotive industry. 
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Nomenclature 
Tabt 300 Ambient temperature (K) 
Papt 101325 Ambient pressure (Pa) 
hair Ambient enthalpy (J/kg) 
Cratr 1005 Heat coefficient at constant pressure (J/ 
(kgK)) 
orf Heat flow rate (J/s) 
Rupstr Upstream specific enthalpy (J/kg) 
R 287 Ideal gas constant (J/(kgK)) 


Aor f Effective cross area 


Crnr(Otnr) Discharge coefficient 

Oct thr Percentage of throttle body that open 
(%) 

y Ratio of specific heats 

Pyatio Pressure ratio 

Prim Pressure ratio limit 

Pupstr Upstream pressure (Pa) 

Paowstr Downstream pressure (Pa) 


Yupstr_air 
Mor f 
Mair 

thr 

Denr 


Qi 


Qi,conv 


Qi cond 
Q, 

Q2,conv 
Q2,cona 


Tintake 

P, intake 
Rintake 
Min intk 
Mout_intk 
Yin_intk 
Vout_intk 
Vin_intk 
Yintake 
Tin _gas 


Tw int 


Text _gas 
Tw ext 


Tmass 
Kint 

Aint cond 
Aint_conv 
Dint_cona 
Kext 


Aext_cond 
Aext_conv 
D ext_cond 
Cq_int 


Cq_ext 


Cpwall 
Myall 

Ven 
ExtnlFlwVe 
i 

Mpuel 

N 

Sinj 

Pwinj 

Ne 
Cc. 


‘yl 

ps 
Tooling 
Ceng 
LHV 
Mpuel 
Texh 


25 
0.003 
0.125 
0.004 
25 


0.003 
0.125 
0.004 


Upstream species air mass fraction 

Gas mass flow through the orifice 

Air mass flow 

Open angle of the throttle (deg) 
Throttle diameter at opening (m) 

Heat flow from the internal gas to a sp 
ecified wall depth (J/s) 

Heat flow convection from the internal 
gas to the internal wall (J/s) 
Conduction heat transfer rate (J/s) 

Heat transfer rate (J/s) 

Convection heat transfer (J/s) 

Heat flow conduction from the external 
middle portion of the wall to the extern 
al wall (J/s) 

Chamber absolute temperature (K) 
Chamber absolute pressure (Pa) 

Specific enthalpy of the intake manifold 
(J/kg) 

Chamber input mass flow (kg/s) 
Chamber output mass flow (kg/s) 
Chamber input heat flow (J/s) 

Chamber output heat flow (J/s) 

Input species mass fraction 

Species mass fraction in the chamber 
Temperature of the gas inside the cham 
ber (K) 

Temperature of the inside wall of the c 
hamber (K) 

External gas temperature (K) 
Temperature of the external wall of the 
chamber (K) 

Temperature of the thermal mass (K) 
Thermal conductivity (W/(mK)) 

Internal conduction area (m?) 

Internal convection area (m?) 

Internal wall thickness (m) 

External wall thermal conductivity (W/(m 
K)) 

External conduction area (m?) 

External convection area (m?) 

External wall thickness (m) 

Internal convection heat transfer coeffici 
ent (W/(m?K)) 

External convection heat transfer coeffici 
ent (W/(m?K)) 

Wall heat capacity (J/(kgK)) 

Thermal mass (kg) 

Chamber volume (m?) 

External flow velocity 


Fuel mass flow (g/s) 

Engine speed (rpm) 

Fuel injector slope (mg/ms) 
Fuel injector pulse-width (ms) 
Number of engine cylinders 
Crankshaft revolutions per power stroke 
Engine cooling temperature (K) 
Heat capacity 

Fuel lower heating value (J/kg) 
Fuel mass flow (kg/s) 

Exhaust temperature (K) 
Engine brake torque (Nm) 
Engine speed (rad/s) 
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Mintake Intake gas mass flow in the chamber (k 
g/s) 
Va Displaced volume 
Leng Engine load 
Ny Engine volumetric efficiency 
Nindicated Indicated load efficiency 
N friction Friction temperature factor 
Tindicated Indicated torque (Nm) 
Tbrake Brake torque (Nm) 
TFriction Friction torque (Nm) 
ARF Air fuel ratio 
Tsta 293.15 Standard temperature (K) 
Tsta 101325 Standard pressure (Pa) 
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